Rush or relax: migration tactics of a nocturnal insectivore in response to ecological barriers

During their annual migration, avian migrants alternate stopover periods, for refuelling, with migratory flight bouts. We hypothesise that European Nightjars (Caprimulgus europaeus) adapt their daily migration tactics in association with biomes. We tracked the autumn migration of 24 European Nightjars, from breeding populations in Mongolia, Belgium and UK, using GPS-loggers and multi-sensor data loggers. We quantified crepuscular and nocturnal migration and foraging probabilities, as well as daily travel speed and flight altitude during active migration in response to biomes. Nightjars adopt a rush tactic, reflected in high daily travel speed, flight altitude and high migration probabilities at dusk and at night, when travelling through ecological barriers. Migration is slower in semi-open, hospitable biomes. This is reflected in high foraging probabilities at dusk, lower daily travel speed and lower migration probabilities at dusk. Our study shows how nightjars switch migration tactics during autumn migration, and suggest nightjars alternate between feeding and short migratory flight bouts within the same night when travelling through suitable habitats. How this may affect individuals’ fuel stores and whether different biomes provide refuelling opportunities en route remains to be investigated, to understand how future land-use change may affect migration patterns and survival probabilities.


Methods
Field methods. We conducted fieldwork in Mongolia (48. . These sites were selected to include individuals from longitudinal extremes of the global distribution range of the species. We captured nightjars in breeding areas using ultra-fine mist nets (Ecotone, 12 × 3 m) and tape lures 24 , marked each individual with a unique alphanumeric ring and fitted a data logger dorsally between the wings 24,25 . In total, we tagged 114 adult individuals, with 1.8 g Pathtrack Ltd. nanoFix or Biotrack Ltd. PinPoint-40 archival GPS-logger (Belgium and the UK) or a 1.2 g SOI-GDL3pam 26 multi-sensor logger (Mongolia and Belgium; hereafter multisensor logger; for a list of deployments see Supplementary Materials Table S1). The multi-sensor loggers contain sensors that recorded air pressure, ambient light intensity, air temperature and acceleration in 5-min intervals and magnetic field in 4-h intervals. The GPS-loggers record longitude, latitude, and altitude at 24-h intervals at 00:00 GMT, when the birds were likely to be in flight (for specific settings see Table S1).

Estimation of migration routes from multi-sensor loggers.
From the recovered multi-sensor loggers, we estimated migration routes by combining flight activity data with light measurements. Multi-sensor loggers measure flight activity as the sum of the absolute differences in acceleration along the z-axis 27 , which allows the detection of flapping flights 26 . As such, we used an automated k-means classification algorithm to classify activity data into three categories: inactivity, low and high activity. This allowed us to define migratory flight bouts as periods (minimum 60 min) of uninterrupted high activity and the delineation of stopovers as periods of one day or longer where migration was interrupted, i.e. where no migratory flight activity was recorded 28 .
We analysed light measurements using a threshold method to provide daily position estimates for the migration period (sensu 29 ) with the R-package PAMLr 28 . Initial position estimates were obtained by identifying sunset and sunrise events from the light-intensity data. We used Hill-Ekstrom calibration, based on light-intensity measurements from the sedentary wintering period, to model the error in defining these sunrise/sunset events caused by shading of the light sensor 30 . Initial position estimates during periods of stopover were grouped together and treated as a single location. Final position estimates were obtained using an Estelle model in SGAT 31 to refine initial location estimates based on Markov chain Monte Carlo (MCMC) simulations and provide a probability distribution around each estimate (two locations per day). In this model we included the following priors: the location where the tags were deployed, the model describing the error in twilight times, a distribution of probable flight speeds (relaxed gamma distribution of shape = 2.2 and rate = 0.08) and a spatial probability mask where stationary positions over water were not possible. We ran this model for 250 iterations to initiate, before tuning the model with the aforementioned priors (three runs with 300 iterations). Finally, to obtain the final position estimates the model was run for 2000 iterations to ensure convergence. Subsequently, we calculated the 97.5% credible intervals around these final position estimates.
Calculation of migration parameters. We only considered migratory flight bouts during autumn migration because most loggers only partially recorded spring migration. Autumn migration was defined as the period between the departure at the breeding areas and the arrival at the wintering areas 25  www.nature.com/scientificreports/ loggers, we excluded data from stopover days. For days containing migratory flight bouts, activity data (inactivity, low activity, high activity, migrating) were categorized into two binary variables: "migrating" (yes/no) and "activity" (yes/no). From previous studies, we know that high and low activity levels are likely associated with foraging behaviour 25 . Hence, "activity" will subsequently be referred to as "foraging activity". For each position estimate, information on the timing of day, night, and twilight (i.e., sunset and sunrise) were extracted using the "classify_DayTime" function from the R-package RchivalTag 32 . We subdivided the daily activity into four time groups: day (from sunrise to sunset), night (from astronomical dusk to astronomical dawn) dusk (from sunset to astronomical dusk) and dawn (from astronomical dawn to sunrise). Air pressure data recorded by the multi-sensor loggers, measured in hectopascals, were converted to estimated flight altitudes above sea level using the atitudeCALC function in the PAMLr R-package 28  where H = flight altitude above sea level (m), T 0 = temperature at sea level (K), L = lapse rate (temperature change per meter increase in altitude, deg/m), P = pressure recorded by the sensor (HPa), P 0 = pressure at sea level (HPa). We used the standard assumptions of P 0 = 1013.25 HPa, T 0 = 288.15 K and L = -0.0065 deg/m 33 . These data on estimated flight altitude above sea level were subsequently used to calculate daily altitude change by subtracting the daily maximum and minimum estimated flight altitude above sea level.
GPS tracks were visually inspected and periods wherein directional movement was interrupted for at least one day were removed, in order to retain only travel segments 34 . From GPS data, we calculated daily travel speed, defined as the distance travelled between consecutive observed locations in 24 h interval periods (km/day). Second, we calculated flight altitude above ground level by subtracting the elevation of the landscape from the recorded altitude. The altitude of the landscape was derived from the CleanTOPO 30 arc second global bathymetric and terrestrial elevation dataset 35 . Grouping data according to biome. To analyse differences in migration tactics in response to biomes, a global classification of natural communities in a particular region based on dominant or major vegetation types and climate 36,37 , we extracted information from the Terrestrial Ecosystems of the World dataset version 2.0 (Dataset published by World Wildlife Fund) for each position estimate and GPS observation. GPS observations were assigned to biomes using the point sampling tool in QGIS version 3.10.10 38 . To account for the spatial inaccuracy of geolocator position estimates 29 , we determined the most abundant biome in the 97.5% credible interval around each position estimate using the zonal statistics tool in QGIS.
The width of the 97.5% credible interval around geolocation position estimates on active migration days was on average 21.23 ± 9.01 degrees in latitude and 9.30 ± 2.14 degrees in longitude. While these credible intervals are large, the most abundant biome represented on average 61.73 ± 15.34% of the surface area. In the East Atlantic flyway, the following biomes were the most abundant in the credible intervals around the different position estimates: "Temperate Broadleaf and Mixed Forests", "Mediterranean Forest, Woodland and Scrub", "Desert and Xeric Shrubland", "(Sub)Tropical Grassland, Savannah and Shrub" and "(Sub)Tropical Moist Broadleaf Forest". In the East Asia/East African flyway these were: "Temperate Grassland, Savannah and Shrub", "Montane Grassland and Savannah", "Desert and Xeric Shrubland" and "(Sub)Tropical Grassland, Savannah and Shrub". Two position estimates from the East Asia/East African Flyway, one categorized as "Boreal Forest" and one as "(Sub)Tropical Moist Broadleaf Forest" were excluded from further analyses since these were the only recorded position estimates in these biomes. In further analyses we distinguished between two categories of biomes, namely ecological barriers, and hospitable biomes (sensu 39 ). The ecological barriers include the "Desert and Xeric scrub" and the "Mediterranean Forest, Woodland and Scrub" biome as well as the "Tropical Moist Broadleaf Forest" biome which is presumed to be a soft ecological barrier for migrating nightjars 34 . Other biomes consist of semi-open habitats and are presumed hospitable biomes.

Statistical analyses.
From daily activity data, we assessed the probability that individuals were migrating in a particular period of the night in response to biome categories. We constructed a binomial model with migratory activity (binary; migrating or not migrating) as the outcome variable and biome category (binary; hospitable or barrier) and period (categorical: dusk, night and dawn) as the predictor. We also included activity of the previous time step to control for temporal autocorrelation. Subsequently we constructed the same model with foraging activity (binary; foraging or not foraging) as outcome variable. Additionally, we investigated how biome category influenced daily altitude change by constructing a similar model with daily altitude change as outcome variable.
From the GPS data, we assessed variation in daily travel speed and flight altitude in response to biome category. We constructed models containing either daily travel speed or flight altitude as the outcome variable and biome category as the predictor. The flight altitude model included a zero-inflation parameter to allow modelling the probability of excess zeros in the conditional part of the model. To investigate the influence of specific biomes on migratory activity, foraging activity, daily altitude change, daily travel speed and flight altitude, we constructed a second set of models where the predictor of biome category was replaced with the variable containing information on the individual biome (Supplementary Materials).
Individual identity was included in all models as a random intercept. Models were fitted using Generalized Linear Mixed Models (GLMMs) with maximum likelihood using the R package glmmTMB version 1.0.2.1 40 . www.nature.com/scientificreports/ Model misspecification problems were checked using the package DHARMa version 0.3.3.0. 41 and when required a zero-inflated model was fitted to model the probability of excess zeros in the conditional part of the model. The significance of the models was established using a likelihood ratio test (Anova with a 'Chisq' test statistic). We then used post-hoc Tukey's comparisons to test the significance of the difference between pairs of interactions between the main predictors using the R package emmeans version 1.5.1. 42 .

Ethics approval and consent to participate
The Mongolian and Belgian research protocols were approved by the Mongolian (Ministry of Environment and Tourism, license numbers: 06/2564 and 06/2862) and Belgian (Agency for Nature and Forest, license numbers: ANB/BL-FF/V18-00086 and ANB/ BL-FF/19-00087-VB) authorities. In the UK, use of GPS tags was approved by the Special Methods Technical Panel and licensed by the British Trust for Ornithology. All protocols were carried out in accordance with the relevant guidelines and regulations.

Results
We recovered 30 loggers (two multi-sensor loggers in Mongolia, four multi-sensor and nine GPS-loggers in Belgium, and 15 GPS-loggers in the UK; all from males; Supplementary Materials S1) constituting a recovery rate of 26.3%. While comparable to other studies (23%; 34 ) this low recovery rate is probably caused by the late deployment of loggers on non-resident individuals with no fidelity to the study site in Belgium (late August 2018) and a low recovery success due to bad weather conditions during a two-week trapping session in Mongolia (July 2019). Of the 30 recovered loggers, six (one multi-sensor logger from Belgium and five GPS-loggers from the UK) did not contain usable data due to a technical failure. Tracked nightjars from Western Europe travel to their wintering grounds in Central Africa along the East Atlantic flyway, with stopovers in Southern Europe, Northern Africa, and the Sahel region (Fig. 1c). Nightjars breeding in Eastern Mongolia migrate to wintering areas in South-Africa and Zimbabwe along the East Asia/East African flyway, with stopovers in Central Asia, the horn of Africa and East Africa (Fig. 1c). The two Mongolian nightjars travelled a minimum migration distance (distance between stopovers) of on average 11696 ± 487 km and spent an average of 100 ± 1.4 days between departure at the breeding site and arrival at the wintering site. Nightjars migrating from Western Europe spent 62 ± 8.1 days on autumn migration, travelling a minimum migration distance of 8633 ± 887 km.
Flight activity data show that nightjars' migration tactics differ between biomes they are traveling through and the period during which they are traveling. There is strong evidence that migration probability during dusk is significantly higher during the crossing of ecological barriers when compared to the category of biomes presumed to be hospitable (Fig. 2, Table 1). Migration probability during the night and dawn periods is not different between biome categories (Fig. 2, Table 1). Our data additionally show evidence that the probability of foraging activity during dusk is higher in hospitable biomes compared to foraging probabilities in ecological barriers (Fig. 2, Table 1). Conversely, nocturnal and dawn foraging probabilities did not differ between biome categories (Fig. 2, Table 1). (for full model results see Supplementary Materials Table S3).
GPS data and barometric pressure recordings show how ecological barriers are crossed at significantly higher speed and altitude. Our GPS data show strong evidence that nightjars cover larger distances per day (Barrier mean = 275 km/day, Hospitable mean = 195 km/day) and fly at higher altitude above the ground (Barrier mean = 1168 m AGL, Hospitable mean = 576 m AGL) when traversing ecological barriers (Fig. 3, Table 2). Similarly, barometric pressure data show strong evidence that daily altitude change is higher during the crossing of ecological barriers (Barrier mean = 1426 m, Hospitable mean = 864 m) compared to hospitable biomes (Fig. 3, Table 2). (For full model results see Supplementary Materials Table S4).

Discussion
During autumn migration, nightjars rush through ecological barriers while migration intensity is relaxed in more hospitable biomes containing semi-open habitats. Rush migration reflects high daily travel speed, high flight altitude and high migration probabilities at dusk. Relaxed migration is characterised by high foraging probability at dusk, lower daily travel speed, lower flight altitude and lower migration probabilities at dusk, suggesting that nightjars may forage and perform short migratory flight bouts within the same night.
A more detailed analysis of migration tactics within specific biomes (Supplementary Materials) shows that Nightjars in our study rush through the Sahara and the Arabian deserts, which are considered ecological barriers for avian migrants in the Palearctic-African migration system 43 . This rushed migration tactic seems to indicate that nightjars attempt to minimize the time spent in these unfavourable areas (sensu optimal migration theory; 6 ). This confirms the findings of a previous study based on a subset of the GPS tracks represented here, where unequal migration speeds where found throughout the flyway 34 . High daily travel speed is also observed in the Central African Rainforest (CAR), supporting the idea that, although interspersed with open areas of human settlements and farmland, the dense tropical rainforest may also form an ecological barrier for Eurasian migratory birds 34,44,45 .
Our data further suggest that nightjars migrate slower in biomes containing semi-open habitats, such as savannahs, grasslands and shrublands. Nightjars' daily travel speed is lower and are interspersed with periods of inactivity and possible foraging behaviour (Fig. 1a,b). It suggests that nightjars may opportunistically feed during the crepuscular hours when they encounter hospitable biomes, which may provide insectivorous birds with relatively large amounts of food 46 . This suggests that nightjars, when encountering semi-open habitats, may travel short distances between sequential sites and accumulate small amounts of extra fuel, a few hours per day. This may correspond to ample observations of nightjars feeding rapidly in key foraging habitats at dusk and dawn during the breeding period 25,47 . At this point it is unclear whether the combination of foraging and migrating  www.nature.com/scientificreports/ reflects a "hop" tactic (i.e. short distance migrations to predictable and reliable food resources 48 ) or a "fly-andforage" tactic (i.e. combine foraging with moving in the direction of travel 19,21,49 ). Nevertheless, both tactics have been observed in combination with a more conventional stopover tactic (i.e. "jump" migration sensu 50,51 ) for other larger-bodied Palaearctic avian migrants. Osprey (Pandion haliaetus) 13,21 and Eleonora's Falcon (Falco eleonorae) 19 , for example, relax their pace across potential foraging habitat in order to combine daily migrating with refuelling, while stopping over for prolonged periods before and after crossing ecological barriers. The acquired fuel reserves en route can therefore serve as an insurance to reduce the time spent at stopover sites before or after crossing ecological barriers 52,53 .
To limit inferences from our data with limited spatial accuracy 30 , which is inherent to the applied geolocation method 54 , we accounted for the error margins around position estimates in assigning environmental characteristics to position estimates and we opted to investigate large scale biomes as opposed to environmental data with a finer spatial resolution. Nevertheless, we acknowledge this method is not without potential biases in evaluating the relationship between environmental characteristics and migration tactics. Future studies investigating how environmental factors influence migration ecology of aerial migrants would therefore benefit from employing a multi-scale approach 55 which combines tracking technologies with fine spatial resolutions such as GPS-tracking 56 .
In accordance with our findings based on data from multi-sensor loggers, GPS-observations in our study indicate nightjars adjust their daily travel speed and flight altitude in response to biomes. GPS-tracked individuals migrating from Western Europe cover approximately 320 km per day across ecological barriers while more suitable areas are crossed at approximately 200 km per day, reflecting how the interplay of topography, wind patterns and food availability influence travel speed of nightjars 57,58 . The observed travel speeds are in line with previous studies on nightjars 34,58,59 and are similar to those of other insectivorous long-distance migrants such as Common Cuckoo (Cuculus canorus) and Common Swift (Apus apus) 60 . Migratory birds are known to adapt not only travel speed but also flight altitude in order to reduce energetic costs and water loss during migration 61 . During the autumn crossing of the Sahara, lower altitudes offer favourable tailwinds 62 and many migratory species, including nightjars, may therefore fly at low altitudes in order to minimize energy costs instead of minimizing water  Table 1.  63 . We observed a mean flight altitude of approximately 1.7 km above the ground when individuals crossed the Sahara or Arabian Desert, similar to other recent studies investigating flight altitude of small nocturnal migrants using barometer logging 26,33,64 . This indicates that when nightjars cross ecological barriers with long flight bouts they climb to higher altitudes, providing better convective cooling, a lack of turbulence and reduced air density 65 , which allows them to achieve higher flight speeds. The lower flight altitude we observe in hospitable biomes may be related to the shorter duration of flight bouts, meaning nightjars were more likely to be at or near ground level during the night when GPS fixes were taken, but also confirm our findings that migrating nightjars are less likely to spend time on active migration when traversing hospitable biomes.

Conclusions
Nightjars adopt variable migrations tactics during active autumn migration. Ecological barriers, where opportunities for foraging are limited, are crossed at a higher pace. This allows individuals to avoid risks and minimize the time they spend in these inhospitable areas. In contrast, where food resources are available, nightjar seemingly spend less time migrating per day while foraging efforts are increased. The ability to opportunistically exploit available resources in between active migratory flight, in our view, could serve to (partially) negate unforeseen poor foraging opportunities en route. While this study provides a first look into how nightjars balance active migration with potential refuelling outside of stationary stopover periods, in future research we will investigate how these different migration tactics during active migration relate to behaviour during prolonged stopovers, allowing investigations on the impact of global land use change on individual migration patterns and tactics.  License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.